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Abstract 

We investigate the collective excitations in spin-one color superconductors. We classify the Nambu-Goldstone modes 
by the pattern of spontaneous symmetry breaking, and then use the Ginzburg-Landau theory to derive their dispersion 
relations. These soft modes play an important role for the low-energy dynamics of the system such as the transport 
phenomena and hence are relevant for late-stage evolution of neutron stars. In the case of the color-spin-locking phase, 
we use a functional technique to obtain the low-energy effective action for the physical Nambu-Goldstone bosons that 
survive after gauging the color symmetry. 
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1. Introduction 

In the phase diagram of quantum chromodynamics (QCD), cold quark matter is expected to deconfine at high 
density and be in a color superconducting state (see for recent reviews). Understanding the properties of matter 
under such extreme conditions is, apart from being an integral part of the quest for the fundamental laws of nature, 
relevant for the astrophysics of compact stellar objects. The density in their cores is presumably high enough to form 
deconfined quark matter. 

Since the energy scale for strong interaction in astrophysical processes mentioned above is of order 100 MeV, one 
can safely neglect the heavy quark flavors and only consider the three lightest ones, up (u), down (d), and strange (s). 
At densities so high that the masses of all light quarks can be neglected, the ground state of three-flavor quark matter 
is known to be the color-flavor-locking (CFL) state. However, as the chemical potential drops to the range interesting 
for astrophysical applications, the strange quark mass is not negligible and starts to play an important role. It reduces 
the Fermi sea of s quarks, leading to an excess electric charge which must in turn be compensated by an imbalance 
between the u and d quarks. As a consequence, the highly symmetric CFL state feels stress and can give way to other 
pairing patterns @]. 

If the stress on the CFL pairing induced by the strange quark mass is too high only the u and d quarks pair, 
which is usually denoted as the 2SC phase. However, as explained above, the energy gain from the 2SC pairing is 
reduced by the requirement of electric charge neutrality. It may then happen that the mismatch between the Fermi 
levels is so large that pairing between quarks of different flavors is completely ruled out. In such a case, one is still 
left with the possibility of pairing quarks of the same flavor. Nevertheless, since the QCD-induced effective quark- 
quark interaction is attractive in the color-antisymmetric channel, this requires the wave function of the Cooper pair 
to be symmetric in Lorentz indices by the Pauli principle, i.e. to carry nonzero spin. While spin one is the simplest 
possibility, the ground state will in general be a mixture of all partial waves with an odd spin |fj. 

There are several situations in which spin-one pairing may occur. Apart from the pairing of quarks of a single 
flavor considered originally in Refs. 17l 4lOll . there is another possibility of pairing of quarks of the same color ll ill . 
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This naturally appears in the 2SC phase where only quarks of two colors are involved in the pairing and the third one 
is left over. We will use this proposal as a warm-up exercise in Sec. [2] However, in the rest of the paper, we will focus 
on the more likely pattern of pairing of quarks of the same flavor. Even then one can distinguish different scenarios. 
When cross-flavor pairing is completely prohibited, all three flavors will pair in some of the spin-one states. On the 
other hand, only s quarks may undergo spin-one pairing as a complement to the 2SC state. We will have in mind 
mostly the first case, since in the latter the presence of the color-asymmetric 2SC state can modify the ground state of 
the spin-one phase 111 211 . 

The classification of the phases of a single-flavor spin-one color superconductor was worked out in Refs. lfl3l [l4ll . 
The physical properties of the so-called inert spin-one phases were investigated in Refs. Il5l4l7ll . Astrophysical 
implications of the presence of a spin-one phase were discussed in Refs. 1 18-25|]. An alternative approach to spin-one 
color superconductors, based on the Schwinger-Dyson equations, was taken in Ref. J26ll . For some recent references 
on the topic see Refs. 1 27 -2^1. 

It was shown in Ref. 11411 that as a consequence of weak interactions, the true ground state of spin-one color 
superconductors is an inhomogeneous state with helical ordering, similar to what happens in non-center-symmetric 
ferromagnets |3(1 31]. In the present paper this additional structure will not, for the sake of simplicity, be considered. 
This may be understood as restricting to length scales much larger than the average distance between quarks, yet much 
smaller than the wavelength of the helix. 

The plan of the pap er is as follows. In Sec. [2] we shall investigate the collective modes in a single-color spin- 
one superconductor [11]. This elucidates some of the peculiarities of spontaneous breaking of spacetime symmetry 
1I32I1 without the technical complications brought by a complex matrix order parameter. In particular, we will show 
an explicit example of a phase in which the number of Nambu-Goldstone (NG) bosons is smaller than the number 
of broken symmetry generators 113 311 (see Ref. 13411 for a recent review including more references). Also, it will be 
demonstrated that due to the fact that the order parameter breaks rotational symmetry, the division of the NG modes 
into multiplets of the unbroken symmetry holds only in the long-wavelength limit. Nonzero momentum of the mode 
breaks the remaining symmetry and results in further splitting of the multiplets. 

Section[3]is the main content of the paper. We analyze the spectrum of collective modes in a spin-one color super- 
conductor, elaborating on the classification of the NG modes suggested in our previous paper [35]. For simplicity, we 
consider quark matter composed of one quark flavor only. This assumption is released in Sec. |4]where we concentrate 
on the color-spin-locking (CSL) phase in electrically neutral three-flavor quark matter. Using a generalization of a 
trick due to Son 1 36], we construct the low-energy effective Lagrangian for the physical NG bosons that remain in the 
spectrum after gauging the color and electromagnetic sectors of the symmetry group. Finally, in Sec. |5]we summarize 
and make conclusions. 



2. Single-color spin-one superconductor 



The pairing pattern considered here was suggested in Ref. [11]. Since it complements the standard 2SC state, it 
involves quarks of the two lightest flavors and one color. They pair in a flavor-singlet state, and the wave function 
must therefore be symmetric in the Lorentz indices. In the simplest case of spin one, the order parameter is a flavor 
singlet and spin triplet, and thus is represented by a complex vector of the SO(3) rotational group. 

We will analyze the excitation spectrum using a Lagrangian which can be regarded as the time-dependent Ginz- 
burg-Landau (GL) theory. We will demand that the Lagrangian has rotational as well as U(l) phase invariance, cor- 
responding to the conservation of particle number, and preserves parity. The most general Lagrangian for a complex 
vector field that has the required symmetries and includes operators of canonical dimension up to four, reads 

«5f = icitf ■ d <p + c 2 <9o0 + • <5 O - a\d^ ■ dtf - a 2 \d ■ <f>\ 2 - b<f> ] ■ <f> - di{tf ■ <pf - d 2 \<f> ■ <f>\ 2 , (1) 

where the spatial vectors are in boldface and their inner product is indicated by a dot; we will use this convention 
throughout the paper. One should keep in mind that by rescaling the field appropriately, one can get rid of one of 
the unknown coefficients, say, c\. The static part of this Lagrangian was investigated in ITlll . so we just summarize 
the results here. First of all, boundedness of the potential from below demands that d\ > and d\ + dz > 0. Once 
b < 0, the scalar field condenses, that is, develops nonzero vacuum expectation value. For d 2 < the ground state 
has the form 0o = v '(0, 0, 1) T , where v 2 = - 2 ( d ^+d 2 ) • ^ e ca ^ ^is the polar phase in analogy with the three-color 



spin-one color superconductor. For d 2 > the ground state can be chosen as 0o = -^=(1, i, 0) r with v 2 = - This is 
analogous to the A-phase. 

In both phases, the SO(3) x U(l) global symmetry is broken to a U(l)' subgroup. In the polar phase, this simply 
corresponds to rotations in the (<p\ , <p 2 ) plane, while in the A-phase, one has to change the overall phase simultaneously 
with the rotation to keep the order parameter unchanged. In the polar phase, all Noether charges of the global SO(3) 
symmetry are zero in the ground state. As a consequence, we expect to find three NG bosons with linear dispersion 
relations at low momentum, associated with the three spontaneously broken generators. On the other hand, in the 
A-state the third component of spin, represented in the 0-space by the matrix = -is gk, h as nonzero density. 



Since /3 belongs to a non-Abelian symmetry group, the NG boson counting will be modified 1331 137H . We expect one 
type-I NG boson with a linear dispersion relation and one type-II NG boson with a quadratic dispersion relation. The 
latter represents a circularly polarized spin wave, very much like in nonrelativistic ferromagnets. We will now see 
how these predictions, based on general properties of spontaneous symmetry breaking in many -body systems [34], 
are verified in an explicit calculation. 

2.1. Polar phase 

The unbroken U(l)' ~ SO(2) subgroup is generated by J3 while the broken generators are J\ .2 and the generator 
of phase transformations. The six-(real-)component complex field <p can therefore be parameterized as 

= <■-•■<■-* J (fa, +H + i X ), (2) 

where 9 and n = (7Ti,n 2 ,Q) T are the NG modes, whereas the remaining degrees of freedom,^ = (~X2,Xi, 0) r and 
H = (0, 0, H) T , are anticipated to excite massive states in the spectrum. The parameterization of the vector ;f is just 
for convenience: to first order in the fields the vector i\ is then "z times 7?' '. Plugging this parameterization into the 
Lagrangian ([]]), expanding up to second order in the fields, and throwing away all total derivatives, we end up with 
the bilinear Lagrangian 

if bilin = -2 Cl {Hd 9 -x ■ d n) + 'Z ) + + (W 2 + (d^xf\ 

-a 2 [(&}H - roUTr) 2 + (d 3 - mt ± x) 2 ] ~ 4(<*i + d 2 )v 2 H 2 + 4d 2 v z x 2 , (3) 

where the planar rotation operator is defined as iot ± tp = d\ipz — d 2 <p\ ■ 

Note that only the H and x modes have mass terms as expected, moreover, the x mass vanishes for d 2 = 0. This 
signals the instability associated with the phase transition from the polar to the A-phase. On the other hand, the H mass 
term can never change sign since d\ + &i > is required by the stability of the potential. All six modes are in general 
mixed by the derivative terms so that a direct diagonalization of the Lagrangian would have to be done numerically. 
However, given that spatial rotations are spontaneously broken to the SO(2) subgroup, we can investigate separately 
excitations along the symmetry -breaking ($3) axis, and in the "transverse" (<p\,(p2) plane. The analysis then simplifies 
considerably and essentially reduces to mixing of two fields at a time. This is discussed in general in | Appendix A| 

Let us first assume that the fields depend just on the time and the third coordinate. The complicated mixing in the 
«2 term then reduces to a mere modification of the phase velocities of the H and 9 fields. Only the singlets H and 9, 
and vectors^ and n now mix. The dispersion relations are found using Eqs. dA.2l > and (lA.3b . In the (H, 9) sector they 
read 

t 1 o t 7 4v 2 (d\ + do)(a\ + a?) n 

E 2 = — [4v 2 (di + d 2 )c 2 + c 2 ] (massive mode), E 2 = ^ — ^-k\ (NG mode), (4) 

c\ 4v 2 (di + d 2 )c 2 + c\ 

to the lowest order in the momentum expansion, while in the (x, n) sector they are 

o 1 9 9 — 4v 2 d 2 ai 9 
E = —r(-4v d 2 c 2 + Cj) (massive modes), E = =ki (NG modes). (5) 

-4v 2 0?2C2 + c 2 

Next we assume that the fields depend just on time and the first two coordinates. The H and 9 modes then decouple 
from others, giving rise to dispersion relations 



1 9 9 9 4v (d\ + d 2 )ai , 

E 2 = -r[4v 2 (di + d 2 )c 2 + c?] (massive mode), E 2 = , v -^-^k 2 (NG mode), (6) 

c\ 4v 2 (di + d 2 )c 2 + c 2 



where the subscript "J." indicates that the mode is transverse. On the other hand, the a 2 term still seemingly mixes 
the components of the ^ and jt vectors. However, redefining them to \ = (x\^Xi.^ T and n - (n 2 , —n\, Q) T , one gets 
r°t±X - d -x an d rot ± 7r = d-fr, while the c\ term is not affected since \ ■ ^ - X' n - The excitations in the (x, n) sector 
then further split into two branches. Modes for which \ an d * are parallel to the momentum behave as longitudinal 
and their dispersion relations are given by 

i 1 t o —4v 2 d*>(a\ + ai) t 

E 2 = —(-4v 2 d 2 c 2 + c\) (massive mode), E 2 = k\ (NG mode). (7) 

Cj -4v 2 d 2 C2 + c\ 

For^ and n modes perpendicular to the momentum, the a 2 term vanishes and the dispersion relations are identical to 
those found in Eq. ((5). 

Note that the masses of the massive modes naturally do not depend on the chosen direction of momentum, as could 
have been expected. Also, it is now obvious that nonzero momentum lifts degeneracy based on the symmetry of the 
ground state. Based on the unbroken SO(2) symmetry we would have expected both^f and n to transform as vectors 
whereas H and as singlets. This is indeed the case for the longitudinal excitations. On the other hand, transverse 
momentum breaks the remaining SO(2) symmetry so that the ni,n 2 NG modes are no longer degenerate. A similar 
remark applies to all other results in this and the following section. 

2.2. A-phase 

The ground state is an eigenvector of ^3, that is, J^q = 0o- This implies that the unbroken subgroup is generated 
by the combination |(1 - J3). The broken generators can be conveniently chosen as J\, J 2 , and |(1 + J3). So <f> can 
be parameterized as 

<f> = e L *" J '(4>o + ^4> +xfa), (8) 

where /' = (Ji,J 2 , ^y 1 ), n = (n\,n 2 ,n^) T , and <f>\ = -^=(1,— /, 0) r is a vector perpendicular to 0o- The real field H 
and the complex field x describe the non-NG modes. When used in the Lagrangian (Q]l, this parameterization yields 
the following bilinear terms 



bilin 



+ y / H = : c i 



f 
a 2 



j(Vi.2) 2 + (d,^) 2 + (d,H) 2 + |5^| 2 



y [(^3^2 - d 2 7T 3 + diX\ + d 2 x 2 + d\H) 2 + (d 3 7Ti - + 82X1 - d\Xi ~ d 2 H) 2 j 
-4d lV 2 H 2 -4d 2 v 2 \x\ 2 , (9) 

where xi,2 ar e the real and imaginary parts of x- Again, thanks to the exponential parameterization the NG fields are 
explicitly eliminated from the static part of the Lagrangian. Also, the^- mass term is seen to vanish at the transition to 
the polar phase, i.e., d 2 = 0. 

To determine the excitation spectrum, we will first assume that the fields depend only on time and the third 
coordinate. The fields then fall into three sectors that do not mix with each other. The field x carries nonzero charge 
of the unbroken U(l)' symmetry and describes a particle-antiparticle pair with masses 



E = ^- (±ci + yjl6v 2 d 2 c 2 + c^ . 



2c 2 

The H and ni, modes give rise, similarly to the polar phase, to one massive and one NG state of type-I, 



(10) 



E 2 = l-(4v 2 dic 2 + c?) (massive mode), E 2 = -^-^±—k\ (NG mode). (11) 

cj 4v 2 d l c 2 +c 2 
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On the contrary, the (n\ , ni) sector, which naively contains two NG bosons, only produces one NG particle of type-II, 
in agreement with the general discussion in 13711 . 



I c l| «l + «2 2 

E- — (massive mode), E — ki (type-II NG mode). (12) 

C2 \Cl\ 

Second, we will investigate the transverse excitations. The masses of the particles excited by Xi>X2 are still given 
by Eq. (fTUl i; the anisotropy brought by nonzero momentum does not appear in this lowest-order term in the dispersion 
relation. For the (H, ^3) fields we obtain the following dispersions, 

1 4v 2 di (ai + f ) 

E 2 = — (4v 2 dic 2 +c?) (massive mode), E 2 = — - -^-k\ (type-II NG mode). (13) 

c 2 Av 2 d\C2 + 

The n\ and 712 modes trivially decouple and give rise to the dispersions 

\C\\ Cl\ j 

E— — (massive mode), E — — k, (type-II NG mode). (14) 
C2 \c\\ 

Before we conclude the section we remark that the calculation of the dispersion relations is complicated by the 
fact that spacetime symmetry is spontaneously broken. However, the basic anticipated features of the NG spectrum 
are preserved. In both phases three generators are spontaneously broken. In the polar phase, they give rise to three 
type-I NG bosons with linear dispersion relation at low momentum. In the A-phase, the ground state carries nonzero 
spin density, therefore the three broken generators produce one type-I NG boson with linear dispersion and one type-II 
NG boson with quadratic dispersion at low momentum. 

Finally, at the transition point between the two phases di = 0, the static part of the Lagrangian has an extended 
SO(6) symmetry under which the polar and A-phase order parameters are degenerate. Five of its generators are broken, 
leaving an SO(5) invariant subgroup. This extended symmetry is also reflected in the NG spectrum. For instance, in 
the polar phase the phase velocities of two of the NG bosons go to zero so that their dispersions are quadratic. This 
is in accordance with the general Nielsen-Chadha counting rule 113 311 . However, note that this extended symmetry is 
explicitly broken by the c\ and (12 terms in the Lagrangian ([TJ. The extra NG bosons are therefore only present in the 
classical theory, they will acquire nonzero masses via radiative corrections 13811 . 



3. Single-flavor spin-one color superconductor 

In this section we will study the spin-one color superconductor that involves pairing of quarks of three colors but 
a single flavor. The diquark condensate or the order parameter A is then a color antitriplet and spin triplet, so it is a 
3x3 complex matrix and transforms as 



U\R 



(15) 



where U = exp(i9 a A a ) e U(3)l = SU(3) C x U(1)b and R = exp(/a,7,) e SO(3)r are transformation matrices. Here 
A a are eight Gell-Mann matrices and /Iq = -Jfl, (Jdjk = — fey* are generators of SO(3)r, a (a = 0, . . . , 8) and a, 
(2 = 1, 2, 3) are rotation angles in U(3)l and SO(3)r group space. There are 18 real parameters in A, among which 12 
parameters are carried by the U(3)l x SO(3)r transformation making a 12-dimensional degenerate vacuum manifold. 
Then A can be parameterized by the remaining 6 real parameters which characterize different vacuum states 111411 . 



Ai 163 —i&2 
-163 A2 i5\ 
162 -id] A3 



(16) 
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3.1. Ginzburg-Landaufree energy and ground states 

The GL analysis is similar to that for superfiuid Helium 3 [39]. Up to fourth order in A and two derivatives, the 
most general U(3)l X SO(3)r and parity invariant Ginzburg-Landau free energy density functional can be written as 



HA] = aiTr(d,A<9jA 1 ') + a 2 (<9 i A fll )((9 ; A* i ) + fcTr(AA t ) 

+t/i[Tr(AA t )] 2 + <f 2 Tr(AA t AA t ) + rf 3 Tr[AA r (AA r ) t ]. 

The time-dependent GL functional, or Lagrangian, is then in general written as 

JS? = iciTrfA^oA] + c 2 Tr[(3 A t )(3 A)] - T[\], 



(17) 



(18) 



see | Appendix B| for more details and the relation between the coefficients c\, c 2 . The ground state is found by min- 
imizing !F[A]. The sign of b determines whether the order parameter is zero or nonzero. Hereafter we will assume 
that b < 0. The phase structure, or orientation in the field space, of Aq depends on d\ and ^3. The magnitude of the 
condensate, v = VTr(AA-f), is given by v = -\J-b/(2d) with d — d\ + /(c/2,^3) (/ is a function of d 2 and d-} that is 
specific to a particular phase) and the vacuum energy is E vac - -b 2 /(4d). The boundedness of the potential from 
below demands that d > 0, which constrains the possible values of d\ for a given phase. The symmetry of the problem 
allows for altogether 8 inequivalent states with different patterns of spontaneous breaking of the continuous symmetry. 
However, only the following four of them occupy a part of the phase diagram [ 14], 
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where a = ^{d 2 + di)/[2(2d 2 + di)] and f3 = ^dtjCldt + c/3). The pattern of spontaneous symmetry breaking de- 
termines the low-energy spectrum of the system, i.e., the NG bosons. While some of the NG bosons are associated 
with the generators of the color SU(3) C group and are thus eventually absorbed in gluons via the Higgs-Anderson 
mechanism, those stemming from spontaneous breaking of baryon number or rotation symmetry remain in the spec- 
trum as physical soft modes. As we will now see, some of the phases exhibit the unusual type-II NG bosons, in 
accordance with general properties of spontaneously broken symmetries in quantum many -body systems l33l l34ll . 
We did not solve the fully coupled equations of motion for the fields depending simultaneously on all coordinates, 
so the dispersion relations shown in the following should be understood as combinations of separate formulas for the 
"longitudinal" and "transverse" excitations. 

3.2. CSL phase 

When d2 + di > and di > c/3, the ground state is the CSL phase whose order parameter is given in Eq. (fT9l . 
Also, d = d\ + ^y 1 . In the CSL phase the spin and color are coupled in the pairing, so the symmetry breaking 

pattern is U(3)l X SO(3)r — > SO(3)v- The generators of the unbroken symmetry SO(3)v are J~^(Aj ® 1 + 1 ® J\), 
yj^(-As ® 1 + 1 ® /2X -\J^(^2 ® 1 + 1 ® /3). There are 9 broken generators leading to 9 NG bosons as follows, 

• Aq®1. Type-I NG singlet, E 2 ~ (a\ + a 2 )k 2 . 

• ^(^7 ® 1 - 1 ® Ji) , Jli-As ® 1 - 1 ® Jz), ^(^2 ® 1 - 1 ® / 3 ). Type-I NG triplet, E 2 ~ (a x + a 2 )k 2 . 

• A a ®l, a = 1,3,4,6,8. Type-I NG 5-plet, £ 2 ~ (a x +a 2 )k 2 . 



The details of the above modes are given in Appendix B. 1 After gauging the color symmetry, the type-I NG 5-plet is 
absorbed by gluons, and we are left with the singlet as the only physical NG boson, stemming from the spontaneous 
breaking of the U(1)b symmetry. 
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3.3. Polar phase 

When dj, < and d 2 + dj < 0, the ground state is the polar phase whose order parameter is shown in Eq. ( fT9] l, and 
d = d\ + d2 + d^. The symmetry breaking pattern is U(3)l X SO(3)r — > U(2)l x SO(2)r. The unbroken symmetry is 
generated by M,2,3 ® 1, ® 1 and 1 ® J$, where P x2 = -^=( V2io + /Is) = diag(l, 1, 0) is the projector onto the first 
two colors. The diquark spin is polarized to one direction. There are 7 broken generators which, however, give rise 
only to 5 NG bosons, organized in the following multiplets, 

• V2!P3 ® 1, where !P 3 = -^(/lo - V2^ 8 ) = diag(0, 0, 1) is the projector onto the third color. Type-I NG singlet, 
E 2 ~ a x k\ + (a x + aj)k\. 

• 1 ® Jj, j = 1,2. Type-I NG doublet, E 2 ~ (a x + a 2 )k 2 ± + a x k\. 

• X a ® 1, a = 4,5, 6, 7. Type-II NG doublet, E 2 ~ a\k]_ + (a x + a 2 fk\. 

The presence of type-II NG bosons is due to nonzero color density of the polar ground state. The details of the above 



modes are given in Appendix B.2 After gauging the color symmetry, the type-I NG singlet and the type-II NG 



doublet are absorbed by gluons, only the type-I NG doublet survives, corresponding to two linearly polarized spin 



3.4. A-phase 

When d-$ > and d 2 < 0, the ground state is the A-phase whose order parameter is shown in Eq. ( [T9| >, and 
d = d\ + d 2 . The symmetry breaking pattern is U(3)l X SO(3)r — > U(2)l x SO(2)y. The residual symmetry is 

generated by A 1^3 ® 1, ^12 ® 1 and -^|(^ > 3 ® 1 - 1 ® J?,). Unlike in the polar phase, the diquark spin is now circularly 
polarized. Among 7 broken generators, there is only one giving rise to a type-I NG mode, 



; 5CP3 ® 1 + 1 ® / 3 ). Type-I NG singlet, E 1 ~ {a x + a 2 )k l ± + a x k\. 
The rest 6 generators produce only 3 type-II NG bosons due to non-zero color and spin density of the A-phase vacuum, 
• X a ® 1, a = 4, 5, 6, 7. Type-II NG doublet, E 2 ~ (a x + a 2 ) 2 k\ + a\k\. 



1 ® ^=(7i ± U 2 ). Type-II NG singlet, E 2 ~ a\k\ + (a x + a 2 ) 2 k\. 



See Appendix B.3 for details on the above modes. After gauging the color symmetry, the type-I NG singlet and the 



type-II NG doublet are absorbed by gluons and leave the type-II NG singlet giving a circular spin wave. 
3.5. e-phase 

When dj > d 2 > 0, the ground state is the s phase, see Eq. (TT~9b for its order parameter; in this case d = 
d x + d 2 ^f - The symmetry breaking pattern is U(3) L x SO(3) R -> U(1) L x SO(2) v . The spin of the second diquark 
color is longitudinally polarized, while that of third color is circularly polarized. Like in the A-phase, the circularly 
polarized spin produces an SO(2)y unbroken symmetry. The unbroken symmetry is generated by y/2'P x ® 1 and 



2 CP 3 ® 1 - 1 ® /3X where f P x = diag(l, 0, 0) is the projector onto the first color. Out of the 10 broken generators 



3 

only two correspond to type-I NG modes: 

'\/2P 2 ® 1, where P 2 = diag(0, 1, 0) is the projector onto the second color. Type-I NG singlet, E 2 ~ a x k\ + (a x 
a 2 )k 2 



, |(!P 3 ® 1 + 1 ® / 3 ). Type-I NG singlet, E 1 ~ (a x + a 2 )k{ + a x k\. 
The remaining 8 generators give rise to 4 type-II NG modes due to non-zero color and spin density of the s vacuum, 
• -^=0li ± iA 2 ) ® 1. Type-II NG singlet, E 2 ~ a\k\ + (a x + a 2 ) 2 k\. 
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-t=0*4,6 ± ^5,7) ® 1. Type-II NG doublet, £ 2 ~ (m + a 2 ) 2 k]_ + a\k\. 
1 ® ± z7 2 ). Type-II NG singlet, £ 2 ~ a\k\ + {a x + ai) 2 k\. 



See Appendix B.4 for details on the above modes. After gauging the color symmetry, only the last type-II NG singlet 
survives, corresponding to a circularly polarized spin wave. 



4. Low energy effective field theory for the CSL phase 

The CSL phase plays in many respects a distinguished role among all spin-one phases investigated in the preceding 
section. First, it is the ground state of one-flavor quark matter in the limit of very high chemical potential. Second, 
it is isotropic and involves democratically all quark colors. Consequently, under some additional assumptions on the 
spin-orbital structure of the order parameter, all quarks are gapped and the low-energy spectrum is solely determined 
by the NG bosons of the spontaneously broken symmetry. Moreover, the rotational symmetry is unbroken whereas 
the NG bosons associated with the color symmetry are absorbed into gluons once this is gauged. Therefore the low- 
energy physics of the CSL phase will be governed by the NG bosons of the global Abelian symmetry. In this section, 
we will employ a technique of Son [36], which was previously used to study the transport properties of the CFL phase 
& 

4.1. Global symmetry and NG bosons 

Based on the argument in the previous paragraph, one would naively conclude that in the one-flavor CSL phase, 
there is exactly one physical NG boson, stemming from the spontaneous breaking of the global U(l) symmetry. 
However, this is only true provided the quarks do not carry any other gauge degrees of freedom apart from the color 
ones. This is certainly a very rough approximation considering that the spin-one phases are only likely to occur in the 
phase diagram in the region where strange quark mass is large and electric charge neutrality effects play an important 
role. 

As already discussed in the Introduction, there are two realistic scenarios for spin-one pairing to occur in three- 
flavor quark matter [1]. (i) The u and d quarks are paired in the 2SC phase and only the s quarks are left and undergo 
the single-flavor spin-one pairing, (ii) Cross-flavor pairing is completely prohibited by large Fermi level mismatch 
and all three quark flavors undergo single-flavor pairing. While in the case (ii) the most favored pairing pattern will be 
CSL, in the case (i) it may not be so. The reason is that the 2SC phase is not color neutral; to compensate for its color 
charge, a color chemical potential must be introduced. This in turn breaks the exact color symmetry of the s quark 
sector. In addition, the mismatch between different colors may favor other pairing patterns such as the polar one [12]. 

For these reasons, we will only consider the scenario (ii) which is theoretically clean; each flavor feels an exact 
color SU(3) symmetry and features the same symmetry breaking pattern. The question then is: how many NG bosons 
are there? As mentioned above, the NG bosons of the spontaneously broken color symmetry will be absorbed into 
gluons. In addition, there is a global Abelian symmetry, G = U(l) i( x U(l) rf x U(l) s , corresponding simply to separate 
conservation of the flavor quark numbers. This will give rise to three NG bosons. However, a subgroup of G, given 
by the electromagnetic U(1)q, is gauged and the associated NG boson will be eaten by the photon, making it massive 
and thus giving rise to the Meissner effect. Therefore, there will be only two physical NG bosons left. 

One should note that in special limits, the flavor symmetry can be actually larger than G. For example, assuming 
that the quark masses satisfy + m u — mj + m s , the flavor symmetry will be G = SU(2)v x U(l) u+d x U(l) s . 
However, the non-Abelian SU(2)v group will be explicitly broken down to U(l)i,, generated by the third component 
of isospin, by the electric charge chemical potential necessary to maintain electric charge neutrality. The remaining 
exact symmetry group is isomorphic to G. The very fact that u and d quarks are actually light does not play a role 
since the pions will presumably still be heavier than the CSL gap so that they will not enter the low-energy effective 
field theory whose validity is limited by the scale of the gap. 
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4.2. General method to construct the effective Lagrangian 

A general method to construct the low-energy effective action for the NG boson of a spontaneously broken U(l) 
symmetry at zero temperature was proposed by Son based on a functional technique lf36ll . Starting with the equation 
of state, P(p), the effective action for the NG field ip to the lowest order in derivatives reads 

r\M, <p\ = j d 4 x p( ^{dw - fi) 2 - (d ilfi ) 2 ). (20) 

Several remarks are in order here. First, this is a fully quantum effective action, that is, it should be used strictly at 
tree level. All loop effects are included in the couplings of the action. Second, after expansion in powers of tp, the 
action contains only terms with the same number of derivatives as is the power of tp. This is the leading-order term 
in the derivative expansion for a given power of ip. Third, the derivation relies heavily on the fact that fi is the only 
source of Lorentz violation in the theory, since then the full dependence on the NG field can be reconstructed from 
the dependence of the pressure on the chemical potential using the fact that the NG field and the chemical potential 
only appear in the effective action in the combination d v tp - 6 v ofi. 

Once we know the equation of state, we plug it into Eq. (1201 and expand in powers of ip to obtain both the dispersion 
relation of the NG boson and its self-interactions. In extremely dense quark matter one can, thanks to asymptotic 
freedom, take as a reasonable starting point the equation of state for a free massless Fermi gas PoG") = N c fi 4 1 (Yin 1 ), 
for a single quark flavor Do, 14011 . The effect of pairing on the equation of state can be neglected since the CSL gap is 
numerically very small. However, one may be interested in corrections due to nonzero quark mass since it is exactly 
the strange quark mass that opens the way to the spin-one phases in the phase diagram. To that end, one needs to 
know the equation of state of a massive Fermi gas, 

Pm(M) Mf wVf , m 4 n + k ¥ fx 4 rn 2 fi 2 4 

— = 12? ' — + log — * 12^ " + 0(m log mX (21) 

where kp = -y'// 2 - m 2 is the Fermi momentum. Substituting this equation of state in the general formula ( f20b . one 
obtains the effective Lagrangian 

i- J^ eff (^) = J-j L 4 - VcW + 6a< 2 (cW) 2 - lAdwf ~ tyto&vtt&V + (d^ip) 2 ] - 
Nc U7T 2 (22) 

- £^ [v 2 - 2/idw + (d 0l p) 2 - (dap) 2 ] . 

The bilinear terms in the Lagrangian imply that the NG boson phase velocity equals 

? 1 2u 2 - 3m 2 1/ m 2 \ ^, A , , ,„„ s 

v =3W^*-A l ~i*r 0{mlii) - (23) 

This can be shown to coincide, to the order displayed, with the hydrodynamic speed of sound in a free gas. 

4.3. Effective Lagrangian for neutral quark matter 

Let us now consider the case of several U(l) symmetries with different chemical potentials. Thus, the sponta- 
neously broken symmetry group G = U(l)i x U(l)2 x ■ • • is associated with the NG fields ip\,<pz, . . . , chemical poten- 
tials fti,fJ.2, ■ • ■ , and the equation of state P{/u\,fj.2, ■ ■ ■)■ Unfortunately, one can easily see that Son's trick does not 
work in this case. The reason is that with more fields there are more independent ways to construct a Lorentz-invariant 
Lagrangian density that reduces to the same function of chemical potentials for uniform fields. As an example, just 
observe that D t 'tpiD fI ip\D v ip2Dyip2 and D fl ipiD v ipiD^tp2D Y ip2 both give n\n\. 

Fortunately, there is a special case where Son's method can still be used. Once the equation of state separates to 

P(MuM2,...) = PM) + P2(M2) + -- - , (24) 

the dependence of the effective action on the NG fields ip^ can again be fully reconstructed using Lorentz invariance. 
This generalization of the method may seem somewhat trivial, since the equation of state d24T > corresponds to separate 
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and noninteracting subsystems. However, they can become entangled by a coupling to an additional field, as we will 
see later. 

Let us address the specific question: what happens when a part of the symmetry group is gauged? For simplicity 
we will assume that there is only one gauge field that couples to a linear combination of generators of U(l)fc, that 
is, to a subgroup of G. This is determined by the charges q k of the fields <p k . In the effective Lagrangian we thus 
have to replace the combinations d v (p k - Syo/Uk with D Y (pk - 6 v ofik, where D^tfk = d^tfk - eqkA^ and e is the gauge 
coupling. To obtain the effective Lagrangian from the equation of state, one in turn has to replace everywhere f3 with 
(Dotpk - Hk) 2 - (Difk) 2 - Note that the effective Lagrangian also contains a term AqJq: a coupling of the gauge field to 
an external background charge density which ensures that the system as a whole remains neutral despite the chemical 



potentials fik This is equivalently expressed by the fact that (Aq) = 0. 

We will now consider a system where the underlying equation of state is well approximated by a noninteracting 
Fermi gas. This is the case of color superconductors at high baryon density since the pairing effects are exponentially 
suppressed and the normal Fermi liquid contribution dominates the pressure. In accordance with Eq. (12211 . the effective 
Lagrangian with the leading finite-mass correction then reads „Sf e ff = Af c (_£?o + -5? i)+_Sf g , where (upon omitting terms 
of zeroth and first order in the fields) 

Jfo =y^2 Z M(A)«fc) 2 - 1$Pm? ~ 4fi k D <p k Dn<p k D fl <Pk + (D ll <f k D u <p k f\ , 

m 2 (25) 

-Sfi = - 2 4^ V D w«) 2 - (Dm) 2 ] , = -4 V v - 

k 

So far we have not fixed the gauge for A M . Note that the covariant derivative d^tfk - eqkA M is invariant under the 
gauge transformation, <p' k — tp k + eq k a, A' fl = A fl + d M a, where a is a gauge parameter. We can change variables ip k to 
9k by (fik = Rkc&c, where Rkc is a real square matrix such that Rki ~ qk up to a common factor (the detailed form of this 
matrix will be specified later). The gauge transformation reads Rff — R6 + eaq or ff — 9 + eaR l q where R,R l are 
matrices and 9', 9, q vectors. We can fix the gauge by choosing the gauge parameter a = - e(R 6 \ hiqk so that 9' x - 0. We 
see that by fixing the gauge we can remove the field 9\, or we can set 9\ - 0. All other modes 62,63,.. . remain in the 
theory as physical fields. The resulting Lagrangian will be rather complicated, let us therefore look explicitly at least 
at the bilinear part of Jz?o 



~£ bilin — / ] 



k {D m f- M * , k {D m f 



(26) 



4tt 2 r \2n 2 
Expanding the square of the covariant derivative one obtains (no summation over k or fi is implied) 

(D^ifk) 2 = (Rkcd^e - eqkA^f = RuRkmd^ed^m - leq^uAyd^t + e 2 q 2 k A fl A fl . (27) 

We can see that the gauge boson acquires a mass term, as expected. Whether this affects the low-energy dynamics 
of the system is a matter of scales. In order to have a clear physical interpretation of the excitation spectrum, it would 
be better to get rid of the mixing term A^d^9(. One could find the dispersion relations even with such mixing, but most 
likely just numerically 1I42I1 . Another, practical aspect is that once we decide to integrate out the massive gauge boson 
to obtain an effective Lagrangian for the NG bosons only, in the absence of the mixing this can be done perturbatively 
and will merely result in a modification of the NG boson interactions. On the other hand, the mixing with the gauge 
boson would necessarily induce corrections to the NG boson dispersion relations. Unfortunately it seems that in the 
most general case of Eq. d2"6*i >. the mixing term cannot be removed by a judicious choice of Ru simultaneously in 
the temporal and spatial parts of Jzfbiiin- However, there are special cases in which the Lagrangian can be further 
simplified. 

(i) Zero masses. In this case, the NG-gauge boson mixing can be removed by choosing Ru so that 2i tAfli.Ru - 

for all I = 2,3, This means that the second and other columns of Ru should be set orthogonal to the vector n 2 k qk, 

that is, not to the linear combination that defines 9\ . It is more convenient to define R k t = [ikRkt since the above 
condition then demands that the ( = 2, 3, . . . columns of R k t be orthogonal to the first one. We can then choose the 
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whole matrix to be orthogonal, upon which the bilinear Lagrangian d26l i becomes 

^biiin = ^ [jJidoOd 2 ~ \idi6tf] + J] e 2 q y k (A A Q - ^A,A,)|. (28) 

l+\ ' k 

The interaction terms are obtained upon expressing tp^ in terms of 6( in Eq. d25b . 

(ii) Equal masses and chemical potentials. In this (rather unphysical) case the mixing term is removed by setting 

Tik QkRke = for all ( = 2, 3, We can thus choose directly the matrix Rkc as orthogonal and the resulting bilinear 

Lagrangian reads 

^bm„ = 2A * 2 4 ~ 2 W2 {X [(cW 2 - vHdiOef] +e 2 Y J ll ( A oA - v 2 A,A,)}, (29) 

t*\ k 

where the phase velocity v is given by Eq. ( l23l . 

All the general formulas above are easily applied to the case of interest, that is, three-flavor quark matter with the 
electric charge neutrality constraint. Then the gauged subgroup is U(1) Q , generated by the electric charge operator, 
Q = (2/3, -1/3, -1/3). While the u and d quark masses can certainly be neglected, the strange quark mass must be 
taken into account, at least for the very reason that without this mass, no mismatch between the Fermi momenta of 
different quark flavors would arise and the ground state would be the CFL phase. 

5. Conclusions 

We analyzed the low-energy physics of spin-one color superconductors in terms of their soft, NG excitations. We 
used the Ginzburg-Landau theory to derive the excitation spectrum. As a warm-up exercise we first analyzed pairing 
of quarks of two flavors and a single color, which may play a role in phenomenology as a complement to the 2SC 
pairing. Already this simple example exhibits a phase with an unusual number of NG bosons. In particular, one of 
the NG bosons in the A-phase is of type-II, i.e., has a quadratic dispersion relation at low momentum, very much 
like the spin waves in ferromagnets. Also, thanks to the fact that a spacetime (rotational) symmetry is spontaneously 
broken we observed that the classification of NG modes into multiplets of unbroken symmetry holds strictly only 
in the long-wavelength limit. Any nonzero momentum of the soft mode further breaks the symmetry and lifts the 
degeneracy based on the symmetry of the ground state itself. 

The main body of the paper is comprised of an investigation of a single-flavor three-color superconductor. This 
is the most likely candidate phase for the ground state of dense quark matter in case that strange-quark mass and 
electric charge neutrality effects disfavor pairing of quarks of different flavors. We have thus completed the analysis 
of the phase diagram started in Ref. II 1 4fl and the classification of the soft modes, sketched in our previous paper 13511 . 
The four phases that appear in the phase diagram possess a plethora of different NG modes of the spontaneously 
broken color, baryon number and rotational symmetry. Those stemming from the color symmetry will eventually be 
absorbed into gluons, making them massive by the Anderson-Higgs mechanism. The other NG bosons will remain 
in the spectrum as physical soft modes. 

Unlike in all the other phases, in the isotropic CSL phase all quarks can be gapped so that the NG bosons are 
the only truly gapless states in the spectrum. This has far-reaching consequences for the low-energy dynamics of the 
CSL phase such as its transport properties. We laid the foundation for a later economical calculation of the transport 
coefficients of the CSL phase such as the shear viscosity by deriving the low-energy effective field theory for its 
physical NG bosons. We used a functional technique [36] applied to the CFL phase before 14011 and adapted it for the 
present case by introducing several independent chemical potentials and gauging a subgroup of the symmetry group, 
corresponding to the electric charge. The actual calculation goes beyond the scope of this paper and will be a subject 
of our future work. 
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Appendix A. Mixing Lagrangians 

Similar to what has been done in 13711 . let us consider the mixing Lagrangian for two fields n, H of the form 

^mixing = \ [(d 7T) 2 - v 2 (dn) 2 \ + I [(d H) 2 - v 2 (dH) 2 ] - Kn 2 H 2 - {Hdo*. (A. 1) 

Here v is the phase velocity; it will be sufficient to assume that it is common to both fields, even though the result 
is easy to generalize to the case with two different phase velocities. The field H possibly has a mass term and £ is 
the mixing parameter which arises from the term in Eq. ([TJ with a single time derivative. The dispersion relations 
following from the Lagrangian (IA. It are 



E = m + 0(k ), massive mode, (A.2) 

2 2 j-4 4 

£l = J T 1 T2 k2 + r 2 ^3 * 4 + 0( * 6) ' NG mode - (A3) 

Note in particular that when the mass term m is zero, the mixing term £ transforms the two expected (type-I) NG modes 
into one massive mode with mass \%\ and one (type-II) NG mode with the quadratic dispersion relation E = v 2 k 2 /\^\. 



Appendix B. Lagrangian and dispersion relations for spin-one color superconductor 

The low-energy effective Lagrangian for the superconductor is in general only constrained by rotational invariance 
as well as the internal symmetry of the system and, in case of QCD, conservation of parity. Its most general form for 
the single-flavor spin-one color superconductor therefore reads 

_S? = ZciTrtA^oA] + c 2 Tr[(d Q A i )(doA)] - aiTr^A^A 1 ") - a 2 (3,A ai )(3 ; Ay 

-fcTr(AA"'") - d\ [Tr(AA"'")] 2 - ^(AA^A 1 ) - J 3 Tr[AA r (AA r ) t ]. (B.l) 

In the limit that the Cooper pairs form tightly bound molecules, the ground state behaves as their Bose-Einstein 
condensate. The coefficients c\ and co are then related. To see this, note that adding a kinetic term to the free 
energy functional T[&\ can be understood as defining a Hamiltonian that governs the dynamics of the molecules, 
J>i? = C2Tr[(<9oA T )(<9oA)] + !F[A]. The many-body description of the system is then accomplished by adding the 
chemical potential /i via ,34? — > Jff - fi^V where jY is the Noether charge corresponding to the U(1)b symmetry 
(baryon number operator). Integrating out the canonical momenta in order to arrive at a Lagrangian formulation of 
the many-body problem [41], one finds terms with one as well as two time derivatives whose coefficients are related 
by c\ - -2/J.C2- Nevertheless, we will keep them as independent parameters. 

Appendix B.l. CSL phase 

With the knowledge of the broken and unbroken symmetry generators, we can write the order parameter A in the 
following form, 

A = exp(/#)exp^i> fl /l fl j(vAcsL + H)exp(iviJi), (B.2) 

where the summations over a and i are in the ranges a = 1,3,4,6,8 and i = 1,2,3. The matrix field H can be 
parameterized by 



H = ht + ]^p a A a +Xih 



i^i + i Xi h- \^ + ^=^ 8 y 6 - ixi 
12 



(B.3) 



where h, ^1,3,4,6,8^1,2,3 are a ll massive fields. We find the following dispersion relations, 



E, = 



Et = 



Ei = 



-2bc 2 + c 1 . (ai + a 2 )(-2bc 2 + 2c?) . 
5 + 7 5 

Cj ~2b c 2 + c \ c 2 

-2bc 2 {d2 — d 3 ) + 3dc 2 2(a\ + a2)\—bc 2 (d2 — di) + 3dc 2 ] , 
— - + ^ — Jr. 



-2bc 2 Xd2 — d 3 ) + 3cfcjC2 
-2bc 2 (d2 + d 3 ) + 3dc 2 2(oi + 02) [-^2(^2 + ^3) + 3<^c?J 2 



3c|<? 



-2bc^{d 2 + di) + 3dc^c 2 



-2b(ai+a 2 ) j2 (fli + gg) c j fc4 

-2/7C2 + C? + (-2Z>C 2 + C 2)3 ' 

-2b(ai + a 2 )(d 2 - d 3 ) 2 21{a x + a 2 ) 2 c\3 i 

—2bc 2 (d 2 — 03) + 3c/Cj 

-2i(ai + a 2 )(^2 + c/3) fc2 + 27 ( fl i + «2) 2 c^ 3 



[-2foc 2 (t/ 2 - d 3 ) + 3^] 3 



-2bc 2 (d 2 + c/ 3 ) + 3c/c^ [-2fcc 2 (^2 + d 3 ) + 3dc]] 3 



(B.4) 



Appendix B.2. Polar phase 

The order parameter can be written in the following form 



c\p(i(n exp| -iK a A a ] (\'Ap„ hir + H I exp( ir,J, .). 



(B.5) 



where the summations over a and i are in the ranges a — 4, 5, 6, 7 and i = 1,2. The matrix field // can be parameterized 
by 



H 



Xn+ifu Xn + ifn 
X21 + upn X22 + Uf22 
ipi ip2 h ) 



(B.6) 



It comprises of two complex doublets of the unbroken SU(2)l, a real vector p of SO(2)r, and a singlet h. They are all 
expected to excite massive modes. The NG mode 6 and the massive mode h are coupled, their eigenmodes are found 
to be 



F z - 



'-2b 

, c 2 

4bc 2 
4bc 2 - 



~2 \ 



2c 2 
lc 2 ) 



1 - 



4bc 2 - c\ t 



^±k 2 + 01 + U2 k 2 
c 2 ± c 2 3 



—k 2 + ai + fl2 / : 2 
c 2 ± c 2 3 



(B.7) 



where ± refers to the orientation of the condensate in the spin space very much like in Sec. [2] The NG mode v and 
the massive mode p are coupled and give the following spectrum, 



.2 \ 
"i 



'ml 
— + 
c 2 2c\ ) 

,2 



1 + 



El = 



2c 2 m- 2 [0.1+0,22 ai,2 



2c2tn 2 + Cj \ C2 



C 2 + 2C2?«2 ' 

kU^k 



Z^+02 a, \ 
\ c 2 c 2 7 



c 2 



(B.8) 



where we have defined m 2 



2 = The NG modes K fl give 



c? / a 1 o 
1+2^ + 



Ol ,2 , a l + a 2 ^2 
C2 3 
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fll + «2 



(B.9) 



The massive modes %ij an d <Pij are coupled to each other and give the spectrum 

i ± 



2 








c 2 




C2 


" 2c\ 





1 + 



4c 2 m 2 



± 1 



c 2 c 2 J / 



Appendix B.3. A-phase 

The order parameter can be written in the following form 

1 



A = expO'0) exp I -iK a X a | (vA A + H) expO'v,/,), 



(B.10) 



where the summations over a and i are in the ranges a = 4, 5, 6, 7 and ;' = 1,2. The matrix // for massive fields can be 
parameterized as 



y?4 - if5 ~i(<P4 - m) Xn + iXn ^ 
<p 6 - itp 7 -i(<p 6 - i(p 7 ) X21 + 'X22 
h + pi + ip2 ih - i(pi + ip2) 



(B.ll) 



The NG fields K4, k$ and the massive fields $4, <ps are coupled (in exactly the same way the NG fields k^, k 7 and the 
massive fields <p-], as they form doublets of the unbroken SU(2) L ) and give the dispersion relations, 

„2 



F 2 

11 K' 



= 4+2 



a\ + ci2 



a -±^lkl + a -±k 2 , 

c 2 c 2 

*1 + ,^ K , 



2 2 

m 2 C \ 

C2 ~ c\ 



11 + 



2c 2 OT2 



+ 1 



1 ± 



\ c 2 + 2c 2 m: 



a\ + ^2 ^2 + gi_^2 
C2 ± C 2 3 



(B.12) 



where we have used m 2 — ~, — 2h(d j di \ m 2 
and give the eigenmodes, 

El 



The NG mode 6 and the massive modes h,p are coupled 



= 1 + 



2bc 2 



a l + fl 2 fc2 + Ol_^.2 
C 2 - 1 C 2 3 



-2* c\ . 

— + 4 +1 1 

Q C 2 



# = 



2bc 2 - c j J \ c 2 



a\ +a 2 2 a i 



/t ~t"" /t^ I , 



C2 



2^ + 27f 



1± J1 + 



(C2in 2 + \bc2) 2 + 4c 2 C2tn 



Ac\ 



1 ± 



4c| 



The NG modes v [i2 give 



y (c 2 m 2 + Abci) 2 + 4c 2 c 2 ffi 2 + 4c| 

4 = 4 + 2(^ + ^1^^, 

C\ \C2 c 2 j / 



^±^ 2 + ^ 2 V 

C2 C 2 



*i 



fll + «2 
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The massive modes Xij 0. j = 1.2) give 



2 








c 2 




± A 


(J 


c 2 


2c 2 





11+^+1 



1 ± 



c\ + Ac2m 2 l 



a\_ k 2 + gj +fl2 /,2 

C2 ± C 2 3 
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Appendix B.4. e-phase 

The order parameter can be written in the following form 



e\p(i6 2 P 2 ) exp(i6 3 P 3) exp | i-K a A a ) (vA E + H) expO'v,/,) 



(B.14) 



where the summations over a and ;' are in the ranges a = 1 , 2, 4, 5, 6, 7 and i = 1,2. The matrix H for massive fields 
can be parameterized as 



H 



<p 4 - i<ps -i(<P4 - m) 
<p 6 - itpj -i((p 6 - iipi) h 2 
K h 3 +pi+ ip 2 ih 3 - i{pi + ip 2 ) J 



(B.15) 



We can define masses 



/;?; 



m 



2b[did 2 + d 2 (d 2 + d 3 )] 2 
d(2d 2 + d 3 ) ' m3 
4 ^lbd\ ^d 2 (d 2 + d 3 ) 2 



Ab{d\ + d 2 )(d 2 + d 3 ) 



d(2d 2 + d 3 ) 
2b(d\ - d\) 



23 



m, 



d(2d 2 + d 3 ) 
„„2 . 



l 45 



d(2d 2 + d 3 ) 



2bd\ 



4bd 2 (d 2 + d 3 ) 



67 



d(2d 2 + d 3 y p d(2d 2 + d 3 ) ' 
The NG modes 9 2 , 8 3 and the massive modes h 2 , h 3 ,pi,p 2 are coupled and give the following eigenmodes, 



mjj + 4m^ 



ml 3 + 4m 2 +4c\/c\ \ c 2 



lll^ 2 + a ' + a2 / : 2 



m 2 ^ + Am 2 



4c 2 



+ 4+1 + 



4c 2 /c 



2 /„2 



m 2 + m 2 3 



ot2 + m,, + 2c, /c 2 



m 2 3 + 4m 2 +4c 2 /c 2 



a \ + fl 2 ^ 2 + ^2 
c 2 ± c 2 3 



a. 2 fli + fl2 
— 1c, + 

c 2 c 2 



4 



2 2 9 

mf + mf, cf 

+ 4 + 

C 2 



4c 2 



1 + 



2c 2 /c 2 



m 2 + m 23 + 2c 2 / C2 , 



c 2 c 2 J / 



2 9 



r. + _1 

2c 2 2c 2 



1± t 1 + 



2c 2 m 2 



1 + 



c 2 + 2c 2 m2 



Qi +fl 2 fc 2 + ^2 
c 2 c 2 3 



For NG bosons k\^ 2 we obtain the energies, 

El 



= 4+2-^1 + 



oi 2 «i + a 2 2 



c 2 



c 2 



fll + «2 
Cl 
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The NG modes £4,5 and the massive modes ^4,5 are coupled (in the same way as K(,j are coupled with ^6,7) and give 
the dispersion relations 



c\ 
a\ + a 2 



2 2 I 

^45 + 5_ 

c 2 2c 2 



^±^i + ^ 2 , 

C2 C 2 / 



Jfci + 



1 ± A 1 + 



2c 2 m 2 5 



+ 1 + 



c\ + 2c 2 m 2 45 



(B.17) 
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The NG modes v\ i give 

E 2 
E 2 
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